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ABSTRACT 


Numerical experiments are conducted with a hydrostatic primitive equation model 
initialized in a baroclinically unstable state to simulate the passage of cold fronts over the 
ocean. The model includes K-theory planetary boundary layer (PBL) parameterization 
with implicitly defined diffusion coefficients. An adiabatic and inviscid simulation 
provided the control for these experiments. The PBL simulations are integrated 1) with 
Zo held constant at 0.4 m and no heat flux; 2) with sea surface temperature (SST) set 
equal to 0i at t = 0 h; 3) with a 5° warmer SST; 4) with diffusion coefficients set equal to 
1. Horizontal resolution is increased to achieve smaller scale fronts in the inviscid and 
ocean simulations. The frontogenetic effects of shear, tilting, convergence, and the PBL 
on isentropic surfaces are evaluated. 

Relative to the inviscid simulation, the PBL simulations produce reduced 
frontogenesis. Surface heat and momentum fluxes combined with turbulent mixing of 
heat promote the development of a deep, well-mixed layer whose depth is dependent on 
the air-sea temperature difference. The rate of frontogenesis is extremely dependent on 
the parameterization of the PBL, specifically surface roughness lengths. Smaller scale 
fronts were produced during the ocean simulations than the PBL land case. Forcing in all 
simulations is due primarily to shearing deformation initially. As the wave grows in 
amplitude, convergence contributes more to frontogenesis than shear. Other terms in the 
frontogenetic equation become important in the PBL simulations. 
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I. INTRODUCTION 


The effect of the boundary layer on the development and modification of fronts as 
they move over various surface conditions is not well understood. Specifically, the 
interaction of fronts with the ocean has been a standing problem to both synoptic and 
dynamic meteorologists. Because the ocean is an infinite source of energy for the 
atmosphere, fronts over the ocean are very sensitive to the temperature differences 
between the air and the sea. For example, when the air is much colder than the water, 
there will be a large heat flux into the atmosphere with strong mixing. A deep mixed 
layer is expected behind a cold front where the air is generally colder than the water. 
This described frontal structure is very similar to that observed by Fleagle et al. (1988). 
They found very large vertical velocities and frontal scales down to a few kilometers. 
When the air is warmer than the water, very stable conditions with much less vertical 
mixing are expected. 

Boundary layer processes have a significant influence on surface fronts because 
they have the largest gradients at the surface in the absence of turbulent diffusion. 
Williams (1974) produced steady state fronts using simple parameterizations of the 
horizontal and vertical turbulent diffusion of momentum and heat. More realistic frontal 
structure was obtained by Keyser and Anthes (1982) who simulated Eady wave 
frontogenesis with a boundary layer representation. The major features in their results 
include a narrow updraft at the top of the boundary layer at the warm edge of the front, a 
stable layer capping the boundary layer to the rear of the frontal zone, and a slightly 
unstable or neutral lapse rate in the boundary layer behind the front. 

Thompson and Williams (1997) used a hydrostatic primitive equation model with 
both K-theory and second-order closure scheme simulations to study maritime 
frontogenesis. In order to identify the boundary layer processes important to maritime 
frontogenesis, moist processes were not included. The results indicated that strong warm 
and cold fronts formed in the adiabatic and inviscid case, but the vertical motion fields 
were comparatively weak. In both the K-theory and second-order closure scheme 
simulations, the boundary layer in the cold air was highly unstable and deep mixed layers 
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formed with a large generation of turbulence. The largest cross-front temperature 
gradients existed in the frontal zone above the mixed layer. These structures were in 
qualitative agreement with observations of maritime cold fronts over the northwest 
Pacific Ocean. However, the simulated fronts did not have the small spatial scales that 
were observed by Fleagle et al. (1988). 

The boundary layer response to changes in sea surface temperature of a frontal 
case over the Atlantic Ocean in FASTEX was predicted by Thompson et al. (2000) using 
the U.S. Navy Coupled Ocean-Atmospheric Mesoscale Prediction System (COAMPS) 
model (Hodur 1997). Initially the front was over colder water and a very shallow mixed 
layer developed in the cold air. The boundary layer was very stable. The front then 
passed over much warmer water and a deep convective boundary layer rapidly 
developed. The increased vertical mixing damped the intense low-level jet on the warm 
side of the front. 

This study investigates the structure of fronts that form over the ocean and the 
boundary layer influence on those fronts passing over the ocean. The basic model is the 
two-dimension primitive equation model employed by Williams et al. (1992) with a first 
order K-theory boundary layer parameterization following Keyser and Anthes (1982). 
Fronts will be forced in this study with a shearing wind field operating on a horizontal 
temperature gradient. Following Thompson and Williams (1997) condensation effects 
will be neglected. Williams et al. (1981) found that condensation did not have a strong 
influence on surface fronts. The secondary goal is to determine frontal structures as a 
function of the air-sea temperature difference. Fleagle et al. (1988) and Chien et al. 
(2001) have observed cold fronts over the northern Pacific with frontal scales as low as 1- 
5 km. These fronts had very little vertical tilt and there was intense mixing in the cold air 
behind the front. Thompson and Williams (1997) were able to simulate this basic frontal 
structure, but could not simulate the small scales because the three dimensional model 
had a grid size of 40 km. 

The simulations are run for both inviscid and planetary boundary layer (PBL) 
modes. The basic model and PBL formulations are developed in Chapter II. 
Frontogenetic forcing terms used to quantify causes of frontogenesis and frontolysis are 
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developed in Chapter III. Numerieal solutions are presented and discussed in Chapter IV, 
and the results are summarized in Chapter V. 
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II. MODEL DESCRIPTION 


This study uses a two-dimensional hydrostatic primitive equation model, whose 
equations are solved numerically using finite differencing on an Arakawa B grid 
(Arakawa and Lamb 1977). Because the vertical scales are much smaller than the 
horizontal resolution, the hydrostatic equation model is used. Finite differences are 
centered in time and space, with an Euler backward time step inserted every fourth 
iteration to control solution separation. Convective adjustment is applied every time step 
by vertically averaging the potential temperature field in areas of static instability. The 
convective adjustment perhaps could be applied less frequently. In the model the time 
dependent quantities are functions of x and z only. The model domain extends 4000 km 
in the east-west (x) direction and 9 km vertically. The horizontal grid spacing varies and 
there are 50 uniformly spaced vertical levels in z. The upper boundary is a rigid lid and 
the horizontal boundaries are periodic. The time step varies depending on horizontal grid 
spacing. 

A. BASIC EQUATIONS 

The equations developed here, following Williams et al. (1992) are adapted to flat 
topography. The vertical axis is normalized with the height of the upper boundary, H, 
and follows 


<7 = 


z 

H 


( 1 ) 


The lower boundary is at o = 0 and the upper boundary is at o = 1. The Boussinesq 
equations (Ogura and Philips 1962) are transformed with equation (1) to give 


av 


av 

+ V«VV + d-— = -V^-ykxV + F 
a<7 


—+ V*V6> + (7—= Q 

dt a<7 


V*V + —= 0 
d<J 


( 2 ) 

(3) 

(4) 
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where 


^ = S.He 

da Oo 


(5) 


0 = t(—) -0o 

\ p y 

</> = [cA +f-c.] 

P 0o 

F and Q are fourth order diffusion terms for momentum and potential temperature plus 
the required boundary layer terms. 

B, INITIAL AND BOUNDARY CONDITIONS 

The frontogenesis simulations are initialized with a small amplitude temperature 
perturbation. The basic state consisting of a 37 m s’' westerly jet at the upper boundary 
decreasing (linearly in z) to 2 m s"' at the surface and an initial temperature field with a 
constant static stability and a small perturbation which is independent of z: 

0{x, z,0) = ^ (z - ^ cos(^ - B) (6) 

dz 2 L 

The V component of the wind is zero to start. The initial temperature perturbation has a 
wavelength equal to the model's horizontal domain, L. Given such a baroclinically 
unstable initial state, rapid frontogenesis may be anticipated. Figure 1 gives the initial u 
and theta fields for all model runs. Constants are contained in Table 1. 

For all but one experiment, the sea surface temperature is set equal to the air 
temperature at the lowest model level at the initial time. It is fixed at this value for the 
duration of the model run. For one run, the sea surface temperature is set to 287 K which 
is approximately 5 degrees warmer than the other simulations. 
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C. BOUNDARY LAYER PARAMETERAZION 
1. Surface Fluxes 

The lower boundary conditions are relatively simple since cases are run 
over the ocean. The parameterization scheme for the surface fluxes is based on Monin- 
Obukhov (1954) similarity theory which gives 


2 

uu* 



2 

vu* 



(7) 

( 8 ) 


w' 0' = —KU * 0 * 


where 


fC 


u* = ■ 


2 2 

U “’“v 


In 


^ z'^ 


y/m 


0 * = 


\Zo) 
K^0l — 0sea^ 


In 


^ z ^ 


y Zoh J 


-y/h 


(9) 


( 10 ) 


( 11 ) 


Here Zo is the surface roughness length, Zoh is the surface roughness length for heat, \|/m is 
the surface layer stability correction function for momentum and \j/h is the surface layer 
stability correction function for heat. The roughness length is calculated using 

=(12) 

exp(-^) 2 

^DN 

Cdn is the drag coefficient for momentum in neutral conditions and K is the von Karmen 
constant. The ratio of the drag coefficients for momentum and heat, R, was estimated by 
Businger et a/. (1971) and is used to calculate Zoh- Initial u* and 0* are calculated using 
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2 2 

Wm ‘*“V 


k{ 0\ — Ose(^ 


In equation (12), the neutral drag coefficient, Cdn, is computed from Vi 
using an algebraic expression (Garrett 1977, 1992) 

Cdn = (0.75 + 0.067Vi)x 10'^ Vi < 20 m s'^ (15) 

Cdn = (0.775 + 0.066Vi)x 10'^ Vi > 20 m s'^ (16) 

where Vi is the wind speed at the 1st level of the model. 

Following Monin-Obovkov (1954) similarity theory, the scaling length, L, 
which characterizes the bouyancy in the surface layer is computed 


where 0 is the average potential temperature at the first level of the model. Given the 
gradient profile functions 

= ~ (IS) 

L u • az 

V|th and t|/m can be determined through integration. The dimensionless similarity gradient 
functions, and d>h, have been deduced from observations. 

In unstable conditions (z/L < 0), the Dyer (1974) and Hicks (1970) 
formulations are used 


1 ^ 
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(21) 


^A = il-\6W 


and \|/h and Vj/mare taken from Paulson (1970) 


,z n (i+k) 

= — - 2a tan(j;) + log 2—^^ 
L Z o 


z J ) 

l^*(-) = 21ogj— 


with y 


--{-kI 


(1988) 


In stable eonditions (z/L > 0), d>m has the form indicated by Hogstrom 


d>„(-) = i + 5- 
L L 


and d>h is determined through the Ellison and Turner (1959) relationship for the ratio —- 


^‘(f) = (l + 4|) 


tj/h and t|/m are taken from Holtslag and De Bruin (1988) 


2 

Z 2 Z 5 -0.35- z 

L 3 L 0.35 L 0.35 


z 2 z 5 -0.35- 2 z . 

L 3 L 0.35 ^ 3 l’ 0.35 
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2. 


K-Theory Boundary Layer Parameterization 

The boundary layer parameterization follows the approach used by Keyser 
and Anthes (1982) which is based on the planetary boundary layer parameterization of 
Blackadar (1979). Turbulent fluxes are represented in K theory by 

- nij 

u'w' = -K.— (28) 

az 


dz 

(29) 

dz 

(30) 


The boundary layer diffusion coefficients and Kh are calculated implicitly from model 
shear and stability fields. The equations are solved in terms of the mixing length,!, 


vertical shear 


dz 


, and the Richardson number, given by 


Ri = 


6o dz 

dv" 

dz 


(31) 


Much debate has been given to what 1 should be. Keyser and Anthes use a constant value 
of 1 = 100 m while most others (Mellor and Yamada 1974; Estoque and Bhumralkar 
1970; Lee and Olfe 1974, Sheih and Moroz 1975; Torrance and Shum 1976) use 


/ = 


A2 


1 + 


/(Z 

J 


(32) 


Louis (1979) used the constant value of X = 100 m, but Thompson and Williams (1997) 
found that X = 150 m for Km and X = 350 m for Kh were more typical for frontogenesis. 
Note that in equation (32) 1 is proportional to height near the surface, but approaches a 
reference value of X at higher elevations. This is physically reasonable since near the 
surface, the size of eddies is constrained by distance to the surface, and 1 should not 
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increase without bounds with elevation. Over the range of slightly negative and positive 
RI, the resulting distribution is approximated by 

K. = Kh = () , Ri > Rie 


= = 1.1 


Ric - Ri 

RU 




dv 


dz 


Ri < Ric 


(33) 


where the eritical Riehardson number is Rh = 0.25 (Blaekadar 1979). Following Keyser 
and Anthes (1982), the eoefficients are broken into a small, eonstant diffusive part Kzo 
and a variable part modeled after equation 33 


K. = K, = A.0 

K. = K, = + 


dv 

f RL-Ri^ 

dz 

1 Rf J 


, Ri > Ric 
, Ri < Ric 


(34) 


where 1 is calculated from equation (32) with X= 150 m for Km and X = 350 m for Kh. 
Finite differencing tends to give greater values of Ri than the mean Ri within the layer 
(Blaekadar 1979), so RIc= 1.0 vice the theoretical value of 0.25 to account for the finite 
grid resolution (following Keyser and Anthes 1982). 

D, FOURTH ORDER DIFFUSION 

In the control run, the model's fourth order diffusion constants are given a value 
of 0.5x10'^ at the surfaee that increases linearly in the vertical to 0.5x10'"^. This linear 
increase maximizes the damping of small-seale numerical noise through the domain 
without giving exeessive surfaee diffusion and it helps to eontrol the reflection of energy 
from the rigid upper boundary. This distribution is not as heavy handed as a traditional 
sponge. 

In the simulations including the planetary boundary layer, the fourth order 
diffusion eoefficients are defined following Keyser and Anthes (1982) as 


K = (Ax)' 


10A10V-r"(Av)' 

2 


du 

ydzj 


+ 




Vdzy 


(35) 
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Horizontal and vertical diffusion coefficients for u, v, and 0 are given equal values within 
the model. 
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III. FRONTOGENETIC FORCING 


The mechanisms that concentrate wind and temperature into a tight zone in a 
frontal region can be explained by examining the individual terms of the frontogenetical 
function, written as 

-= Fadv + Fcon + Ftilt + Fsh + FpBL ( 37 ) 

dxdt 

where 


d 

fdo^ 

d 

fdo^ 

'""dx 

[dxj 

"^dz 

[dxj 


Fcon — 


Ftilt — — 


Fsh = — 


du do 
dx dx 

dw do 

dx dz 

dv do 
dx dy 


(38) 

(39) 

(40) 

(41) 


dx dz 


dx 




d d 

H- 

dx dz 


Kh — 

dx 


|-( 0 ) 

ox 


(42) 


Fadv (38) represents the advection term. The convergence term (39) quantifies the 
effect of u on isentropes in reducing or concentrating the temperature gradient. The 
tilting term (40) quantifies the effect of a gradient of vertical motion tilting the isentropes 
into a vertical plane. Shear deformation (41) quantifies the effect of the v wind on the 
north-south horizontal gradient in the isentropes. Fpbl (42) represents forcing due to 
planetary boundary layer diffusion, 4th order diffusion, surface momentum flux and 
surface heat flux. 
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IV. NUMERICAL RESULTS 


The results from the numerieal simulations are discussed in the following 
sections. The emphasis of this study is on the comparison of the inviscid case results to 
those of the planetary boundary layer simulations which were run with the same initial 
conditions. Table 2, summarizes the simulations conducted with 10 km horizontal 
resolution and At = 30 s. Simulations A and C are run at 5 km horizontal resolution with a 
At = 10 s. 


Simulation 

A 

B 

c 

Cl 

C2 

Momentum 

flux 

No 

Yes 

Yes 

Yes 

Yes 

Zo 

- 

Zo= 0.4 m 

Calculated 

Calculated 

Calculated 

Heat Flux 

No 

No 

Yes 

Yes 

Yes 

SST 

- 

- 

01-282 K 

287 K 

01-282 K 

Zoh 

- 

- 

Calculated 

Calculated 

Calculated 

4th order 

diffusion 

0.5x10'^ 

inereasing 

linearly 

vertical to 

0.5x10'' 

Eq (36) 

Eq(36) 

Eq(36) 

Eq(36) 

PEL 

diffusion 

No 

Eq (35) 

Eq(35) 

Eq(35) 

II 

II 

s o 


Table 2. Numerical Simulations 

A. ADIABATIC AND INVISCID SIMULATION (A) 

The model is run first with the surface flux and boundary layer parameterizations 
excluded to serve as a control case. The case is nearly inviscid in the sense that fourth 
order horizontal diffusion is always present. The initial u and 0 fields are shown in 
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Figure 1 as functions of x and z. The cross frontal velocity field in Figure la (and all 
subsequent figures) includes the mean background flow. 

Due to the baroclinically unstable initial state, rapid frontogenesis may be 
anticipated. Figure 2 contains the frontal solutions after 108 hours of integration. Only 
the lower portion of the atmosphere is shown to highlight the near-surface effects and 
only a 1500 km horizontal extent surrounding the front is shown to closely examine the 
area around the front. The front has intensified and all fields have been advected 
downstream with the background current. The vertical motion field, Figure 2c, shows 
ascending motion all ahead of the front with weak descending motion in the cold air 
behind the front giving a thermally direct circulation about the front. Also the v field 
shows large cyclonic vorticity just east of the maximum temperature gradient (Figure 2b). 


a. Uatt = 108 b. Vat 108 



Figure 2. Inviscid frontal solutions at t = 108 h; (a) cross-front flow, m s'^; (b) along-front 
flow, m s’'; (c) vertical velocity, m s'^; (d) potential temperature, deg K. 
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Rapid frontogenesis continues since the north-south temperature gradient does not 
change. Figure 3 contains the frontal solutions after 120 hours of integration. The front 
has continued to intensify. The v field indicates increased horizontal shear and there is an 
increased horizontal temperature gradient (Figure 3b-c). The front is strongest at the 
surface as indicated by the increased horizontal temperature gradient in that region. The 
southerly wind ahead of the front is bringing in warm air while the northerly wind behind 
the front is bringing down cold air from the north (Figure 3b,d). 


a. U att= 120 b. Vat 120 



1650 1900 2150 2400 2650 2900 3150 1650 1900 2150 2400 2650 2900 3150 


km km 

Figure 3. Inviscid frontal solutions at t = 120 h: (a) cross-front flow, m s'^; (b) along-front 
flow, m s’'; (c) vertical velocity, m s'^; (d) potential temperature, deg K. 


Figure 4 is a comparison of the frontogenetic forcing terms Fsh, Ftiu and Fcon at t = 
108 h and t = 120 h with different scales. Values of the frontogenetic equation indicate 
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that, while the wave is small amplitude, the main forcing is due to Fgh- Opposing this 
frontogenesis is the Ftiu term. Because this term goes to zero at the surface the front is 
strongest at the surface and weaker above. Convergence in the u field has continued to 
become stronger. As the wave grows in amplitude Fcon will contribute more to 
frontogenesis than Fgh- Also note that the maximum values of Fgh and Fcon are to the east 
of the maximum temperature gradient region. This causes the front to move to the east 
relative to the mean flow, which allows it to keep its upwind vertical tilt. 


a. Shear. (Inviscid) at t = 108 
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Figure 4. Frontogenetic forcing terms for inviscid simulation; (a) Fsh at t = 108 h; (b) Fsh at 
t = 120 h; (c) Ftiit at t = 108; (d) Fat at t = 120 h; (e) Fcon at t = 108 h; (f) Fcon at t = 120 h; 
the potential temperature lines (dashed) are included for reference. 


To quantify the scale of the frontal zone, the parameter D is defined as 
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(36) 


D = 


|A^| 

w~ 


dx 


max 


where A0 is the maximum horizontal temperature variation at that level. Beeause D is a 
measure of the temperature gradient and no frietion is eonsidered, it is expected to 
oscillate when no frontogenesis occurs and to collapse when a front forms. The frontal 
scale for both the frictionless and PBL land cases is presented in Figures 5 and 6. Figure 
5 gives a comparison at 90 m and Figure 6 gives a comparison at 1.35 km. At t = 168 h, 
the D = 34.7 km at 90 m and D = 55.5 at 1.35 km. This is due to the fact that the tilting 
term is frontolytic, but goes to zero at the surface. Since the frontolytic effect decreases 
towards the lower boundary the frontal scale is smaller at 90 m. 


Frontal scale atz = 90 m 



Figure 5. Frontal scale parameter at 90 m; dotted black line with o: case A (inviscid); solid 
green line: case B (PBL over land). 
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Frontal scale at z = 1.35 km 



Figure 6. Frontal scale parameter at 1.35 km; dotted black line with o; case A (inviscid); 
solid green line; case B (PBL over land). 

Figure 7 shows a time series of the maximum along-front flow in the lowest 16 
levels of the model (2.79 km) for the present case (A), the K-theory boundary layer land 
case (B) discussed in section 3.2 below, and the K-theory boundary layer ocean case (C) 
described in section 3.3 below. Note the rapid intensification beginning near 72 hours. 
Figure 4 contains the frontal solutions after 120 hours of integration when the front is 
extremely strong. The vertical circulation has increased slightly. The ascent ahead of the 
front at t = 108 h is only 0.027 m s"' (Figure 2c), but by t = 120 it has increased to 0.044 
m s'^ (Figure 3c). At this point the along-front velocity, as shown in Figure 7, is 
becoming unrealistic. Because the north-south temperature gradient, dO/dy, does not 
change, this front will continue to grow. 
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Maximum Along-Front Velocity 



Figure 7. Maximum along-front velocity from the surface to 2.79 km; dotted black line with 
o: case A (inviscid); dashed blue line: case C (PBL over ocean); solid green line: case B 
(PBL over land). 


These results clearly show that in the absence of boundary layer processes, a 
growing baroclinic wave will generate a strong front at the surface. However this 
simulation is unrealistic because along frontal wind velocities exceed 60 m s'^ by t = 126 
h (Figure 7). In addition, the circulation and vertical motion about the front is too weak. 
During the Storm Transfer and Response Experiment (STREX) in November 1980 
vertical motion fields for a cold front were observed by aircraft to be 6 m s"' in the 
updraft and 2 m s'^ in the downdraft (Bond and Eleagle 1985). 

To determine whether a smaller scale front could be achieved, simulation A was 
run at Ax = 5 km. Eigure 8 shows a time series of frontal scale calculated with (36) for 
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both horizontal resolutions. After t = 126 h, frontal seale for the Ax = 5 km solution 
begins to deviate from the Ax = 10 km solution and beeome smaller. By t = 180 h, D = 
11.1 km for Ax = 5 km and at t = 168 h, D = 23.1 km for Ax = 10 km. Clearly smaller 
seale fronts ean be aehieved using a higher resolution model. Similar numerieal results 
were obtained by Williams (1967). 
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Figure 8. Frontal seale at z = 90 m; black dotted with o: Ax = 10 km; black dotted with □: 
Ax = 5 km. 

B, LAND SIMULATION 

For this experiment, the initial conditions are identical to those in the adiabatic 
and inviscid case. A typical value for woodland area of 0.4 m (Garrett 1980; Fichtl and 
McVehil 1970) is assigned to Zq. Since sensible heat flux is small over land except for 
during daylight hours, sensible heat flux is neglected at the lower boundary in this case to 
idealize land without a diurnal cycle. 
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Figure 9 compares the potential temperature fields at t = 6 h for cases A and B. 
By t = 6 h, a well-mixed layer is already developed in the PBL land case (Figure 9b). 
The potential temperature gradient associated with the front is much different in this 
simulation than in the inviscid case. The well-mixed layer continues to develop and 
deepen in the PBL land case through t = 126 h (not shown). 

a. Theta (Inviscid) at t = 6 b. Theta (PBL Land) at t = 6 




Figure 9. Potential temperature fields at t = 6 h: (a) inviscid case; (b) PBL land case. 

Figure 10 presents frontal solutions for the PBL land case at t = 108 h. The 
frontal zones depicted in the potential temperature fields at t = 108 h exhibit no tilt in the 
vertical from the surface to approximately 1 km (Figure lOd). The well-mixed areas on 
both sides of the front prevent tilt. This is in contrast to the inviscid case in which the 
frontal zone has pronounced tilt (Figure 2d). The location of the maximum along-front 
wind speed is no longer at the surface, but at approximately 1 km at t = 108 h (Figure 
10b) with lower wind speeds near the surface. In addition the overall maximum wind 
speed is lower (Figure 7). The lower wind speed near the surface can be attributed to 
surface drag. This case is similar to the cases treated by Keyser and Anthes (1982) and 
Peng et al. (2001). 

The effect of the addition of surface momentum flux is clear from Figures 5 
through 10. Surface momentum flux produces a well-mixed surface layer, reduces the 
along-front velocity and reduces the growth of the baroclinic wave since the frontal scale 
is larger in the PBL land case. In addition at the surface, it takes longer for the 
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frontogenesis to begin and for the PBL front to develop. The PBL frontal scale is smaller 
above the surface as seen by the lower value of D. 


a. Uatt=108 b. Vat 108 



450 700 950 1200 1450 1700 1950 450 700 950 1200 1450 1 700 1950 



Figure 10. PBL frontal solutions at t = 108 h: (a) cross-front flow, m s'*; (b) along-front 
flow, m s'*; (c) vertical velocity, m s'*; (d) potential temperature, deg K. 

The baroclinic wave amplitude is smaller in this case and the potential 
temperature gradients associated with the front are different from the inviscid simulation. 
There are several reasons for this. Looking at the frontogenetic terms, Fpbl becomes 
more important in this simulation since surface momentum flux and PBL diffusion have 
been added. The effect of Fpbl is frontolytic. Figure 11 displays Fsh for the inviscid and 
PBL land simulations at t = 72 h with different scales. This time period was chosen 
because in the inviscid simulation after t = 72 h the front intensifies rapidly. Fsh has the 
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same order of magnitude in each simulation, but for the PBL case the maximum 
frontogenetic effect is not at the surface and is about one half the magnitude of the 
inviscid case. The overall frontogenetic effect is reduced further in the PBL case because 
PBL diffusion is firontolytic. Figures 5 and 6 indicate that the frontal scale for the PBL 
case was smaller at 1.35 km than at 90 m. This is due to the frontogenetic effect of Fgh 
going to zero at the surface since it is proportional to v. 


a. Shear. (Inviscid) at t = 72 xIO ” ^ Shear. (PBL Land) at t = 72 



Figure 11. Fgh at t = 72 h: (a) inviscid simulation; (b) PBL land simulation. 

C. OCEAN SIMULATIONS (C, Cl, C2) 

Three ocean simulations are conducted. In cases C and C2, the sea surface 
temperature is set equal to the air temperature at the lowest model level at the initial time. 
It is fixed at this value for the duration of the model run. For case Cl, the sea surface 
temperature is set to 287 K which is approximately 5 degrees warmer than for the other 
simulations. In case C2, virtually all PBL diffusion is removed by setting Km = Kh = 1.0 
at all grid points except at the lowest level. 

Figure 12 shows a time series of the maximum along-front wind for all ocean 
simulations and the inviscid simulation. As expected in the case with the warmer sea 
surface temperature, the frontal scale is not as small as in the simulation run with the 
colder water due to decreased mixing. Figure 13 illustrates a time series of maximum Km 
for simulations B, C and CL Until t = 80 h. Km values are largest for the land simulation 
indicating the strongest mixing. Km values are larger for the warm water simulation than 
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the cooler water simulation up until t = 80 h. Case C2 with small planetary boundary 
layer diffusion more closely approaches the inviscid case, but because of the addition of 
the surface momentum and heat fluxes it is not as strong. 


Maximum Along-Front Velocity 



Figure 12. Maximum along-front velocity from the surface to 2.79 km: dotted black with o: 
inviscid simulation (case A); dashed blue with x: ocean simulation with small PBL 
diffusion (case C2); dashed blue: ocean simulation with 282K SST (case C); dash-dot 
red: ocean simulation with 287K SST (case Cl). 
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K values 



Figure 13. Km values: green solid: PBL land simulation (case B); dashed blue: ocean 

simulation with 282K SST (case C); dash-dot red: ocean simulation with 287K SST (case 
Cl). 


Figure 14 shows the time evolution for frontal scale as defined by (36) for the 
ocean simulations and the inviscid simulation. Although the along-front velocity in all 
ocean simulations is less than in the inviscid simulation, from t = 72 h until t = 114 h 
(Figure 14) the frontal scale is smaller in the ocean simulations. 
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Frontal Scale at z = 90 m 



Figure 14. Frontal scale at z = 0.9 km; dotted black with o; inviscid simulation (case A); 
dashed blue with x: ocean simulation with no PBL diffusion (case C2); dashed blue: 
ocean simulation with 282K SST (case C); dash-dot red: ocean simulation with 287K 
SST (cased). 

In the ocean simulations it takes longer for the well-mixed layer to develop than 
in the PBL land simulation. At t = 12 h, cases C and C2 are just starting to form a well- 
mixed layer while case Cl already has (not shown). The larger upward heat flux 
associated with the warmer seawater in case Cl is causing more turbulence and stronger 
mixing in the surface layer than in the cooler SST cases. Remember that in the PBL land 
simulation, a well-mixed layer was formed by t = 6 h (Figure 9), but as expected in the 
ocean simulations it takes longer. The delay in the formation of this well-mixed layer can 
be attributed to the dependence of surface heat and momentum fluxes on the surface 
roughness lengths and in these cases Zo (.001 - .002 m) is significantly smaller than the 
0.4 m used in the PBL land simulation. 
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1. Effect of Planetary Boundary Layer Diffusion and Surface Fluxes 

By comparing potential temperature fields for cases C and C2, the timing for 
when boundary layer diffusion (35) begins to have an influence ean be determined. 
Potential temperature fields for these cases are identical up until t = 66 h (not shown) and 
frontal scale measured by D is equal through t = 54 h (Figure 14). By examining eddy 
viscosities values, Km and Kh, as defined by (35), it beeomes apparent that boundary layer 
diffusion only affeets case C after t = 54 h when Km and Kh exeeed 1. The greater values 
of Km and Kh after t = 54 h are consistent with larger values of u*. Figure 15 shows the 
frontal solutions at t = 108 h for cases C and C2. Adding the K-theory boundary layer 
parameterization results in a more realistic 0 pattern. The potential temperature gradient 
aeross the front is larger in C compared to C2 (Figure 15g-h). This potential temperature 
gradient remains well defined in ease C despite boundary layer friction. These 
parameterizations also result in more realistic along-front velocity fields. While C2 is 
more realistie than the inviseid simulation, it still continues to grow and produces 
unrealistic along-front velocities that exeeed 60 m s"' after 5.5 days. The addition of 
boundary layer diffusion in case C produces a reasonable along-front wind veloeity of 36 
m s'^ at day 6. 

Figure 13 indieated that higher heat flux in case Cl produced larger Km values up 
until t = 80 h. After t = 80 h, ease C Km values are larger than ease Cl 

Comparison of Figure 15 to the corresponding frontal solution for the inviseid 
simulation (Figure 2) clearly shows the influence of the boundary layer diffusion, surface 
momentum flux and heat flux. There is clear evidence of the frontal structure in the 
temperature fields for the oeean simulations. The strongly stratified cold air behind the 
front in the inviseid simulation is replaced in these simulations by a well-mixed layer 
approximately 1 km deep for case C and 0.5 km deep for case C2. Strong surface fluxes 
in the unstable air behind the front in case C promote entrainment and a greater 
deepening of the boundary layer than case C2. The frontal zones depicted in the potential 
temperature fields (Figure 15g-h) exhibit no tilt in the vertical over this well-mixed layer. 
This is in contrast to the inviseid simulation, but it agrees with strong fronts observed 
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over the oeean (Bond and Fleagle 1985; Fleagle et al. 1988; Chien et al. 2001). The laek 
of tilt is due to well-mixed boundary layers on both sides of the front. 


a. U (Ocean SST - 282 K) at t = 108 b. U (Ocean SST - 282 K no PBL diffusion) at t = 108 
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c. V (Ocean SST - 282 K) at 108 d. V (Ocean SST - 282 K no PBL diffusion) at 108 
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e. W (Ocean SST - 282 K) at t = 108 Hours f. W (Ocean SST - 282 K no PBL diffusion) at t = 108 Hours 
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Figure 15. Frontal solutions for cases C and C2 for t = 108 h: (a) U (m s"') for ocean 

simulation SST 282 K (case C); (b) U (m s"') for ocean simulation SST 282 K no 
PBL diffusion (case C2); (c) V (m s'^) for case C; (d) V (m s'^) for case C2; (e) W (m s'^) 
for case C; (f) W (m s’') for case C2; (g) Theta (K) for case C; (h) Theta (K) for case C2. 
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2, Frontal Scale 

Figure 16 shows a time series of frontal scale calculated with (36) for both 
horizontal resolutions for case C. After t = 126 h, frontal scale for the Ax = 5 km begins 
to deviate from the Ax = 10 km and become smaller. By t = 168 h, D = 54.5 km for Ax = 
5 km and at t = 168 h, D = 76.2 km for Ax = 10 km. Clearly smaller scale fronts can be 
achieved using a higher resolution model. During the PBL land simulation with Zo= 0.4 
m, the smallest D reached was D = 136.6 km at t = 186 h. The simulations conducted 
produced smaller scale fronts over the ocean than over land. 


Frontal scale at z = 90 m 
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Figure 16. Frontal scale at z = 90 m; blue dashed: Ax = 10 km; blue dashed with A: Ax = 5 
km. 


3, Frontogenetic Forcing in the Ocean Simulations 

As in the inviscid simulation, values of the frontogenetic equation indicate that 
the main forcing in the ocean simulations while the wave is small amplitude is due to Fgh. 

Opposing this frontogenesis is the Ftin term. Convergence in the u field becomes stronger 
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as the wave amplifies. As the wave grows in amplitude, Fcon will contribute more to 
frontogenesis than Fsh. Figure 17 shows Fcon normalized by the absolute value of Fsh for 
the inviscid simulation and the ocean simulation, C. When Fcon is normalized by the 
absolute value of Fsh exceeds 1, Fcon is contributing more to frontogenesis than Fsh. Fcon 
contributes to frontogenesis more than Fjh only after 5.5 days for both cases. For the 
inviscid simulation, Fcon grows significantly faster than in the ocean simulation in which 
frictional effects decrease the wind speed at low levels and increase low-level cross- 
isobaric flow. The frontogenetic effect of Fcon in the inviscid simulation exceeds Fsh 
approximately 6 hours prior than in the ocean simulation. 

25 

20 

15 

10 

5 

1 

0 

Figure 17. Fcon normalized by the absolute value of Fsh: solid with *: inviscid simulation 
(Case A); dashed black: ocean simulation with 282K SST (Case C). 

Other terms in the frontogenetic equation become important in the ocean 
simulations since surface momentum flux, surface heat flux and PBL diffusion have been 
added. Because all simulations are thermally direct circulations, the tilting term is 
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predominantly frontolytic east of the front beeause of the distribution of vertieal motion. 
As the vertical motion fields become more compact and the maximum in vertical motion 
moves closer to the lower boundary the Ftin gets a dipole structure as seen in Figure 18. 
Larger gradients in vertical motion produce larger magnitude in the tilting term. The 
largest gradient in vertical motion is within the frontal zone where the vertical motion 
changes sign. In Figure 18, this area causes a reversal in direction of the gradient in 
vertical motion and results in a frontogenetic contribution from the tilting term. In the 
ocean simulations, after t = 96 h an area of weak sinking develops east of the frontal 
zone. During the ocean simulations, the tilting term increases 3 orders of magnitude in 
the frontal region by t = 144 h. While the tilting term does this in the inviscid simulation, 
by t = 144 h Ftin is still 3 orders of magnitude smaller than the shear term. For the ocean 
simulation by t = 144 h, Ftiu is only one order of magnitude smaller than Fsh. 



Figure 18. Ftiitat t = 102 h for case C; dotted black: W (m s’'). 

4, Effect of Sea Surface Temperature 

In an attempt to determine frontal structures as a function of the air-sea 

temperature, the simulation Cl was conducted with an initial sea surface temperature 
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(SST) of 287 K which is approximately 5 K higher than the SST in simulation C. 
Beeause the air in this simulation is colder than the water, a deeper mixed layer is 
expected behind the cold front then was seen in simulation C. Due to the increased 
turbulent mixing in case Cl less frontogenesis is expected than for ease C. 

Figure 19 shows the frontal solutions at t = 108 h for oases C and Cl. As 
expected the warmer sea surfaoe temperature is causing increased mixing resulting in a 
deeper mix layer. The well-mixed layer is deeper in ease Cl through t = 126 h (not 
shown). After t = 126 h, the well-mixed layer appears to be approximately the same 
height, but for case C it is cooler. The potential temperature gradient aoross the front is 
larger in the case C. The along-front wind velocity is less, and the oonvergenoe at the 
surfaoe is less resulting in less ascending motion for ease Cl. With the warmer sea 
surfaoe temperature cooling near the surfaoe is reduced so the oold pool behind the front 
is warmer than in case C. Beeause the cold pool is cooler in case C, the temperature 
gradient within the inversion at the top of the boundary layer is stronger. 

Maximum along-front wind velooity is less in case Cl than ease C at all times 
(Figure 12) this is due to increase vertical mixing as indicated by higher Km values until t 
= 80 h (Figure 13). The frontal scale for case Cl is larger than case C at after t = 80 h, 
but by t = 168 h both simulations have the same D value. The additional heat flux in ease 
Cl initially slows frontogenesis, but eventually both fronts acheive the same small scale. 
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a. U (Ocean SST ~ 282 K) at t = 108 


b, U (Ocean SST = 287 K) at t = 108 
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e. W (Ocean SST - 282 K) att = 108 Hours 



g. Theta (Ocean SST ~ 282 K) at t = 108 




f. W (Ocean SST = 287 K) at t = 108 Hours 
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Figure 19. Frontal solutions for cases C and Cl for t = 108 h: (a) U (m s"') for ocean 

simulation SST 282 K (case C); (b) U (m s"') for ocean simulation SST = 287 K (case 
Cl); (c) V (m s'^) for case C; (d) V (m s’') for case Cl; (e) W (m s’') for case C; (f) W (m 
s'^) for case Cl; (g) Theta (K) for case C; (h) Theta (K) for case Cl. 
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V. CONCLUSIONS 


The numerical simulations conducted in this study indicate that the PBL plays an 
important role in maritime frontogenesis. The effects in the inviscid simulation are 
consistent with those found in other studies. Specifically, it produced intense fronts with 
strong convergence and vorticity. Frontogenetic forcing was due primarily to shearing 
deformation. Using higher horizontal resolution, smaller scale fronts are achieved. 

Relative to the frictionless results, the PBL simulations produced more realistic 
frontogenesis. Surface heat and momentum fluxes combined with turbulent mixing of 
heat promoted the development of a deep, well-mixed layer. In the PBL simulations the 
fronts do not tilt in the vertical over the lowest portion of the atmosphere. This structure 
is consistent with observations of fronts over the northeast Pacific Ocean (Bond and 
Fleagle 1985). 

Frontogenesis is extremely dependent on the parameterization of the planetary 
boundary layer. Specifically, surface roughness lengths. During the land simulation, 
surface roughness length is fixed at 0.4 m and a well-mixed surface layer forms 
immediately. Frontogenesis is extremely slow and weakened during the land simulation. 
This agrees with cases treated by Keyser and Anthes (1982) and Peng et al. (2001). 
During the ocean simulations where surface roughness lengths are calculated and range 
from 0.001-0.002 m, the well-mixed layer takes longer to form. The addition of surface 
momentum flux clearly produces a well-mixed surface layer, reduces the along-front 
velocity and reduces the growth of the baroclinic wave. It takes longer for the 
frontogenesis to begin and for the PBL front to develop. 

The addition of surface fluxes is not enough to properly simulate frontogenesis. 
Planetary boundary layer diffusion is also necessary. In the inviscid case rapid 
frontogenesis continues because the north-south temperature gradient does not change. 
In the simulation with small PBL diffusion, surface drag reduces the near surface wind 
and the convergence in the boundary layer, but the front still continues to grow and 
produces unrealistic along-front velocities by 5.5 days. 
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Air-sea temperature differenees determine the depth of the well-mixed layer. In 
both oeean simulations as the eold air is brought in from the north, the air beeomes eolder 
than the water. This temperature differenee is inereased for simulation Cl whieh has a 5° 
higher water temperature. Turbulenee and mixing are inereased with the larger upward 
heat flux assoeiated with the warmer seawater in simulation Cl as seen in the higher 
planetary boundary layer diffusion values. 

Frontogenetie foreing in all simulations is due primarily to shearing deformation 
initially. As the wave grows in amplitude, eonvergenee contributes more to frontogenesis 
than shear. Other terms in the frontogenetie equation become important in the PBL 
simulations since surface momentum flux, surface heat flux and PBL diffusion have been 
added. 

Smaller scale fronts were achieved in the ocean simulations than in the land 
simulation agreeing with observational studies (Bond and Fleagle 1985; Fleagle et al. 
1988; Chien et al. 2001). Increased horizontal resolution produced even smaller scale 
fronts. 

In the future, the addition of other surface layer elements including moisture 
should be investigated and higher resolution simulations should be conducted to try and 
achieve smaller scale frontal features. Condensation enhances surface convergence 
locally by warming the column through the release of the latent heat. The impact of this 
associated latent heat release on frontogenesis can be important. Latent heat release, due 
to the lifting of moist unstable air in the lower levels of the troposphere, appeared to be 
critical to frontogenesis. Depending on the level of convection and corresponding latent 
heat release, the vertical motion profiles will be altered and therefore the tilting term of 
the frontogenetie equation will be changed. In addition to areas of increased ascent 
within the clouds there will be areas of descent between clouds. The location of 
maximum cyclonic vorticity is dependent on the vertical motion distribution and 
therefore the distribution of convective heating. Condensation will also redistributing 
surface air upwards so that comparatively high 0e air is located along the tropopause. 
Because the extent and level of condensation can vary significantly the over all effect on 
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the vertical motion field and heating make it difficult to hypothesize on whether the 
addition of moisture would be frontogenetic or firontolytic. 
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